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Abstract. A new class of solvable A^-body problems is identified. They describe N unit- 
mass point particles whose time-evolution, generally taking place in the complex plane, 
is characterized by Newtonian equations of motion "of goldfish type" (acceleration equal 
force, with specific velocity-dependent one-body and two-body forces) featuring several ar- 
bitrary coupling constants. The corresponding initial-value problems are solved by finding 
the eigenvalues of a time-dependent N x N matrix U(t) explicitly defined in terms of the 
initial positions and velocities of the N particles. Some of these models are asymptotically 
isochronous, i.e. in the remote future they become completely periodic with a period T 
independent of the initial data (up to exponentially vanishing corrections). Alternative for- 
i— i mulations of these models, obtained by changing the dependent variables from the N zeros 

^ of a monic polynomial of degree N to its N coefficients, are also exhibited. 
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\q 1 Introduction 

q 

^) In this paper a new class of solvable iV-body problems is identified. They describe an arbitrary 

number N of unit-mass point particles whose time-evolution, generally taking place in the comp- 
lex plane, is characterized by Newtonian equations of motion "of goldfish type" (acceleration 
equal force, with specific velocity-dependent one-body and two-body forces, see below) featuring 
several arbitrary coupling constants. The solvable character of these models is demonstrated by 
the possibility to ascertain their time evolution by purely algebraic operations. In particular it 
is shown below that their initial-value problems are solved by finding the eigenvalues of a time- 
dependent N x N matrix U (t) explicitly defined in terms of the initial positions and velocities 
of the N particles. Some of these models are asymptotically isochronous, i.e. in the remote 
future they become completely periodic with a period T independent of the initial data (up 
to exponentially vanishing corrections). Alternative formulations of these models, obtained by 
changing the dependent variables from the N zeros of a monic polynomial of degree A to its iV 
coefficients, are also exhibited. 

The main idea to obtain these results is to identify the A coordinates z n (t) characterizing 
the positions of the particles of the A-body problem with the A eigenvalues of an A x A matrix 
U(t), itself evolving according to a solvable (or integrable) matrix ODE. This technique was in- 
vented long ago by Olshanetsky and Perelomov [1-4] and has been much exploited subsequently, 
identifying thereby many solvable (or integrable) many-body problems: see for instance [5] (in 
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particular its Section 2.1.3.2 entitled "The technique of solution of Olshanetsky and Perelomov 
(OP)") and [6] (in particular its Section 4.2.2 entitled "Goldfishing"), and the more recent pa- 
pers [7-23]. The present paper provides new results of this kind, by taking as point of departure 
a solvable ODE characterizing the time-evolution of the N x N matrix U(t) which is different 
or more general than those previously employed to this end. These results are treated in the fol- 
lowing section, including its subsections and Appendices A and B where all the equations solved 
in this paper are listed (hence, the reader wishing to get an immediate glance at them may 
immediately jump to these appendices). The last section, entitled "Outlook", tersely outlines 
further developments whose treatment is postponed to subsequent papers. 

2 Solvable iV-body problems 

In this section we describe two classes of solvable A r -body problems. The models of the first 
class are not new, hence their treatment is not elaborated beyond their identification; several 
models of the second class are new, hence they are fully dealt with. 

In Subsection 2.1 we introduce a system of two matrix first-order ordinary differential equa- 
tions (ODEs) characterizing the time-evolution of the two N x N matrices U = U(t) and 
V = V(t), and we indicate how the corresponding initial- value problem can be explicitly solved. 

In Subsection 2.2 we show how - via the introduction of two appropriate ansdtze - the N 
eigenvalues z n (t) of the matrix U(t) can be identified with the N coordinates of N unit-mass 
point particles whose time-evolution, generally taking place in the complex plane, is characterized 
by Newtonian equations of motion ( "acceleration equal force" , with nonlinear one-body and two- 
body forces). These models are thereby shown to be solvable. The first ansatz yields various 
models whose solvability was already known; hence their treatment is limited to the derivation 
of their equations of motion. The second ansatz yields several new models - as well as several 
previously known models - characterized by equations of motion with specific velocity-dependent 
one-body and two-body forces "of goldfish type" featuring several arbitrary coupling constants. 
They are all listed in Appendix A. The alternative class of many-body models obtained by 
changing the dependent variables from the N zeros of a monic polynomial of degree N to its N 
coefficients is discussed in Subsection 2.3; the corresponding equations of motion are listed in 
Appendix B. 

As it is well known (see for instance [5], in particular Chapter 4, entitled "Solvable and/or 
integrable many-body problems in the plane, obtained by complexif ication" ) , models such as 
those described below, which describe the motions of N points in the complex z-plane, can 
be reformulated as iV-body models describing the motion in a plane of N point-particles, the 
positions of which are identified by real 2- vectors in that plane. While we leave the elaboration 
of this connection to the interested reader, we feel justified by it to refer to our findings as 
describing "physical" (if not necessarily "realistic") many-body problems in the plane. 

2.1 A solvable system of two matrix ODEs 

In this subsection we discuss the following system of two matrix ODEs, satisfied by the two 
N x N matrices U and V: 

U = aU + f3U 2 +-/V + rj(UV + VU), V = pV. (2.1) 

Notation. Upper-case Latin letters generally denote N x N matrices (unless otherwise indi- 
cated), with the (scalar!) N being throughout an arbitrary positive integer. Here and below the 
matrices are time-dependent (unless otherwise indicated), in particular U = U(t), V = V(t). 
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Lower case letters generally denote scalars. The 5 scalar parameters a, (3, 7, 77, p are time- 
independent but a priori arbitrary, until specific restrictions on their values are explicitly men- 
tioned. Superimposed dots always indicate time differentiations. 

Remark 2.1. The factor (3 multiplying the term U 2 in the right-hand side of the first of the 
two matrix ODEs (2.1) could be eliminated (i.e., replaced by unity) by rescaling [/(and by 
correspondingly replacing 7 with an adjusted value, say 7); then the factor r\ multiplying the 
term {UV + VU), or the (adjusted) factor 7 multiplying the term V, could also be eliminated 
(i.e., replaced by unity) by rescaling V. It is preferable not to do so in order to keep open the 
possibility to set one or the other of these parameters, n or 7, to zero (but we will never set 
both of them to zero, to avoid decoupling the time evolution of U(t) from that of V(t)). Note 
moreover that the introduction of a constant scalar term (implicitly multiplied by the N x N 
unit matrix /) in the right-hand side of the first equation would amount - up to a redefinition 
of other parameters - to adding to the matrix U the N x N unit matrix / multiplied by a new 
(constant) parameter, implying just a constant shift of all the eigenvalues of the matrix U, 
a trivial change not worth pursuing; while the introduction of two different parameters in front 
of the two terms UV and VU could be eliminated by adding, in the right-hand side of the first 
of the two matrix ODEs (2.1), the commutator [U, V] of U and V times a convenient parameter, 
without any effect on the eigenvalues of U. If some of the parameters in the matrix ODEs 
(2.1) vanish this model may reduce to one of those that have been previously treated, see [6] 
(in particular its Section 4.2.2 entitled "Goldf ishing" ) and the more recent papers [7-23]; we do 
not discard these models in the following, since we consider in any case worthwhile to exploit 
the unified treatment provided by the present approach based on (2.1). As for the possibility 
of replacing the right-hand side of the second matrix ODE (2.1) with a more general function 
of the matrix V, it is a possibility whose exploration is postponed to a future investigation (see 
Section 3). 

Solution of the system of two matrix ODEs (2.1). Clearly the solution of the second of the 
two matrix ODEs (2.1) reads 

V(t) = V exp(pt), V = V(0). (2.2) 

To solve the first of the two matrix ODEs (2.1) it is convenient to set 

U = -^V-(f3F)- 1 F, (2.3) 

obtaining thereby for the N x N matrix F = F(t) the following second-order linear matrix 
ODE: 

F — aF + FW = 0, (2.4a) 

W = {-acq + /3 7 + r}p)V -r/V 2 . (2.4b) 

It is then plain via (2.2) that the general solution of this matrix ODE reads 

F(t) = F+f+(t; V Q ) + F_/_(t; V ), (2.5) 

where F± are two arbitrary constant N x N matrices and the two scalar functions f±(t;v) 
(which of course become N x N matrices when the scalar v is replaced by the N x N matrix Vq, 
see (2.5)) are two independent solutions of the scalar second-order linear ODE 



/ — af + [(— cxr\ + /?7 + rjp)v exp(pt) — rj 2 v 2 exp(2pi)] / = 0. 



(2.6) 



4 



F. Calogero and G. Yi 



And it is easily seen that two independent solutions of this ODE are given by the following 
formulas: 



f+(t;v) = exp 
/_ (t; v) = exp 



nv 

exp(pt) 

P 

— exp(pt) 
P 



^(-^l-^eMpt)), (2.7a) 
V 2r)p p p ) 

eM*m(^^-,l + ^eMpt))- (2.7b) 
V 2r)p p p J 



Here the (scalar!) function &(a; c; z) is the confluent hypergeometric function (see, for in- 
stance, [24]). 

Note that the formula (2.3) with (2.5) implies the following explicit solution of the initial-value 
problem for the matrix U(t): 

U(t) = -i [ V V exp(pi) + [f+(t; V ) + C Vo)]" 1 [/+(t; V ) + C/-(i; V )] } , (2.8a) 

with the time-independent matrix C defined in terms of the initial data U(0) and Vq = V(0) 
(see (2.2)) as follows: 

C = -{/+(0; [0*7(0) + ryF ] + /+(0; Vb)}{/-(0; F o )[/3*7(0) + 77V0] + /_(0; V )}-\ (2.8b) 

Remark 2.2. The fact that the expression (2.3) (with (2.2) and (2.4)) of the matrix U(t) satis- 
fies (2.1), and likewise that the functions f±(t;v) (see (2.7)) satisfy the ODE (2.6), can be easily 
verified. Of course these formulas are valid for generic values of the parameters that appear in 
them, excluding special cases - such as vanishing values of p or 77 (see (2.7)) or of j3 (see (2.8a)), 
or values of a and p such that a/p is an integer (in which case the two functions f± (t;v) do 
not provide two independent solutions of the ODE (2.6): see Section 6.7 of [24]) - which are 
clearly problematic (although the formulas may in these cases be reinterpreted via appropriate 
limiting procedures). We ignore hereafter these issues, even when listing below solvable iV-body 
models a few of which belong to these problematic cases. Indeed in this paper we mainly limit 
our consideration to the identification of solvable iV-body problems; these findings open the 
way to obtaining a rather detailed understanding of their actual behaviors (see for instance the 
following Remark 2.3), but such analyses exceed the scope of this paper: they should be done 
on a case-by-case basis, and shall perhaps be postponed to the moment when some of these 
iV-body models evoke a specific, theoretical or applicative, interest. 

Remark 2.3. If the parameter p is purely imaginary and the parameter a is real and negative, 
27ri 

p=— , a<0 (2.9a) 

(of course with T real and nonvanishing, and i the imaginary unit, i 2 = —1), then clearly the 
matrix U(t) (see (2.8) with (2.7)) is asymptotically isochronous (i.e., asymptotically periodic 
with the period T independent of the initial data): indeed in this case, as t — > +00, 

U(t) = U+(t) + 0[exp(at)] (2.9b) 

with U+(t) given by the formula (2.8) with C = 0, 

U+(t) = -- {rjVoexp(pt) + [f + (t; ^o)] _1 [/+(*; ^0)]} , (2.9c) 

hence (see (2.9a) and (2.7a)) it is periodic with period T, 

U+(t + T) = U+(t). (2.9d) 

This observation entails of course the property of asymptotic isochrony for the solvable N- 
body models identified below featuring parameters a and p which satisfy the conditions (2.9a). 
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2.2 Identification of solvable many-body models 

The starting point is to introduce the N eigenvalues z n {t) of the matrix U(t), and the corre- 
sponding diagonalizing matrix R(t), by setting 

U(t) = R(t)Z(t)[R{t)}-\ Z(t) = diag[z„(t)]. (2.10a) 

Here and hereafter indices such as n, m, £, k run over the integers from 1 to N (unless otherwise 
indicated). 

Likewise we set 

V(t) = R(t)Y(t)[R{t))-\ Y nn {t) = y n (t), (2.10b) 

and we introduce the N x N matrix M(t) by setting 

M(t) = [#(*)] -^(i), M nn (t) = iJ, n (t). (2.10c) 

Remark 2.4. The diagonalizing matrix R(t) is defined up to right-multiplication by an ar- 
bitrary diagonal N x N matrix D(t), since the replacement of R(t) by R(t) = R(t)D(t) does 

not affect (2.10a). But it changes M(t) (see (2.10c)) into M(t) = [Rfr)]' 1 R(t) implying the 
following change of its diagonal elements: fJ, n (t) fi n {t) = ^n{t) + d n (t) / d n (t) , where the N 
quantities d n {t) are the, a priori arbitrary, elements of the diagonal matrix D(t). Hence we 
retain the privilege to assign at our convenience (see below) the diagonal elements /u n (i) of the 
matrix M(t). 

It is now easily seen that via these assignments (2.10) the two matrix ODEs (2.1) get rewritten 
as follows: 

Z + [M,Z]=aZ + /3Z 2 + -fY + 7](ZY + YZ), Y + [M, Y] = pY. (2.11) 

Here and hereafter the notation [A, B] denotes the commutator of the two matrices A and B: 
[A, B] = AB — BA. 

Let us now look, componentwise, at the diagonal and off-diagonal elements of these two 
matrix ODEs. 

The diagonal part of the first of the two ODEs (2.11) reads 

z n = az n + f3zl + -yy n + 2r]z n y n , (2.12a) 
implying 

Vn = (2.12b) 

7 + 2 V z n 

The off-diagonal part of the first of the two ODEs (2.11) reads 

- (z n ~ z m ) M nm = [7 + rj(z n + z m )]Y nm , n^m, (2.13a) 
implying 

7 + 77 (z n + z m ) 



Y nm , n^m. (2.13b) 



Zn Z m 

The diagonal part of the second of the two ODEs (2.11) reads 

N 

y n = pyn+ Yl (YneM en - M ne Y en ) , (2.14a) 
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implying, via (2.13b), 

N 

Vn= PUn + 2 { Y n(Xln 



7 + v(z n + z £ ) 



zi 



(2.14b) 



Hence, by time-differentiation of (2.12a) we see via (2.12b) and (2.14b) that the coordina- 
tes z n (t) satisfy the following system of ./V equations of motion of Newtonian type ( "acceleration 
equal force", with one-body and two-body forces): 



-apz n - Ppz n + (a + p)z n + 2f3z n z n + 

N 



2rjz n (i 



az„ 



P4) 



7 + 2r l z r. 



+ 2( 7 + 2nz n ) \ Y mYin 



7 + V( z n + Zj) 
Zn ~ Z£ 



(2.15) 



But in these equations of motion the role of "two-body coupling constants" is played by the 
quantities Y n gYi n (with £ ^ n) which are in fact time- dependent. Indeed the time evolution of 
the off-diagonal elements Y nm of the matrix Y is determined by the off-diagonal part of the 
second of the two ODEs (2.11) which componentwise read 



N 



Y nm = pY nm + (Y nk M km - M nk Y h 



nk^kra) > 



n 7^ 771, 



(2.16a) 



k=l 



yielding, via (2.10b), (2.10c), (2.12b), (2.13b) and a bit of algebra, 



Y 



Y n 



Z n — Z m Tj (z n OLZ n /3z n ) 7] (z m OLZ m /?z m ) 

a + p + /3{z n + z m ) h ' 



N 



(I, 



+ Mm + 



Zn Z m 



y 

J- 71: 



7 + 2r\z n 
2n + (7 + 2rjzi) 



h — 

Zl Zn 



7 + 2r\z 
1 



n ^ m. 



(2.16b) 



So, in order that (2.15) become the TV" Newtonian equations of motion of a genuine iV-body 
problem one must either provide some "physical interpretation" for the quantities Y n £ (with 
£ ^ n) - possibly in terms of internal degrees of freedom, an alternative we do not pursue in 
this paper - or find a way to "get rid" of these quantities, i.e. find a way to express them - if at 
all possible - via the N coordinates z m = z m (t) and possibly also the N velocities z m = z m (t). 
Previous experience [5-23] suggest that two types of ansdtze are the appropriate starting points 
to try and achieve this goal. 



2.2.1 First ansatz 

The first ansatz reads 



y 



{OnQn 



,1/2 



•"n 



n 7^ 771, 



(2.17) 



where we reserve the option to make a convenient assignment for the N functions g n of the 
coordinate z n , g n = g n {z n ) = 9n[z n (t)]- 

The insertion of this ansatz in (2.16b) yields, after a bit of trivial algebra, 

9lL_^±P_Bz - V ^ n ~ ^ - ?® + ((n -> m)) 
2g n 2 PZn 1 + 2nz n +^ n ^ m )) 
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g n (7 + 2r]z m ) - g m (l + 2rjz n ) 



N 



Zn Z m ) 

7 + 2 VZn 



(z n - zgf 



E {* 



7 + 2rjz m 
(Zm - Zif 



n ^ m. 



(2.18) 



Here and throughout the notation +((n — > m)) indicates the addition of whatever comes before 
it, with the index n replaced by m. 

We now take advantage of the freedom (see Remark 2.4) to assign the diagonal elements p n 
of the matrix M by setting 



N 



7 + 2r?z„ 



{Z n - Zi) 2 _ 

and we moreover make the assignment 
9n = 9(7 + 2r]Z n ), 



(2.19) 



(2.20) 



with g an arbitrary constant (i.e., g = 0). Thereby the system of N(N — 1) equations (2.18) 
gets reduced to the following, much simpler, system of only N algebraic equations: 



a + P a , yjazn+jh%) _ n 
„ P z n + _ — u, 

2 7 + 2r]z n 



(2.21) 



which amounts to the following 3 equations (recall that we exclude the uninteresting possibility 
that 7 and r] both vanish): 



Pv = o, 

$7 + VP = °> 
7 (a + p) = 0. 

This boils down to the following 3 possibilities: 

a = /3 = p = 0, 



or 



or 



/3 = 7 = p = 0, 



P = V = 0, P 



-a. 



(2.22a) 
(2.22b) 
(2.22c) 



(2.23a) 



(2.23b) 



(2.23c) 



The corresponding solvable iV-body models - obtained by inserting (2.17) with (2.20) and with 
these assignments of the parameters in the Newtonian equations of motion (2.15) - read, in the 
first 2 of these 3 cases - after conveniently setting 



7 + 2rjz n (t) = exp[2<„(i)] 
with c an arbitrary nonvanishing constant - as follows: 



Cn 



g 2 rf 



\ d N 

\JL £ s inh- 2 [c(Cn-0)]; 



(2.24a) 



(2.24b) 
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while in the third of these 3 cases they read 

N 



2_,4 



a z n + g 7 



d_ 

dz r , 



^ (Z n - Zi) 



(2.25) 



But these are well-known solvable iV-body problems, see for instance [5]. Hence we conclude 
that from the matrix system (2.1) no new solvable many-body models are obtained via the 
ansatz (2.17). 



2.2.2 Second ansatz 

We proceed then to consider a second ansatz, reading 

Y n m — \_9n9m {^n fn) (Zm ~t~ fm)\ > Tl ^ J71, 



(2.26) 



with g n and f n functions of the coordinate z n that we reserve to assign later: g n = g n (z n ) = 
g n [z n (t)], fn = fn(z n ) = f n [z n (t)]- Its insertion in (2.16b) yields, again after a bit of trivial 
algebra and now the assignment fi n = (again justified by Remark 2.4), the following system 
of N(N - 1) second-order ODEs: 



1 / Z n + fn . 9n 

9 _l_ f a ~ P 



j3z v 



az„ 



(3z n 



N 



+ V(z n + fn)9n ~ ^2 + ft)9t 

Z n [gn{l + 2r\Z n ) - 1] + fn9n(l + 2r]Z n ) 



To reduce this system we clearly must now set 



7 + 2f]Z n 
7 + rj(z n + 



+ ((n — > m)) = 0, n ^ m. 



1 



9n 



7 + 2r]z n ' 

fn = f (0) + f W Zn + f 2) zl 



(2.27) 

(2.28a) 
(2.28b) 



where /^°', f^\ f^ 2 ' are 3 constant parameters that we reserve to assign later. Thereby the 
system of N(N — 1) ODEs (2.27) gets transformed into the following system of (only) N ODEs: 



A' 



Zn 2(z n -|- fn) ^ ^ 



a + p + 2/3z n + 



(_Z£ + fe)h + V(Zn + Zi)] 

(z n - zi){^ + 2r]zi) 

2r] (z n - az n - (3z n - f n ) 
7 + 2rjz n 



+ fn) + {f {1) + 2f (2 h n )f n . (2.29) 



Note that to write this system in more compact form we employed a mixed notation, using 
sometimes the functions f n or fg instead of their explicit expressions, see (2.28b). 

To complete the task of ascertaining for which values of the (so far arbitrary) 8 constant 
parameters a, (3, 7, r/, p, f(°\ f^\ f^ the system of N(N — 1) ODEs (2.29) is satisfied we 
must utilize the equations of motion (2.15) which, via the ansatz (2.26) with (2.28), now read 



N 



2(in + /n) 

£=l,£^n 



(M + fd\l + VjZn + Z £ )] 
(Z n ~ 2<)(7 + 2l]Z£) 
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2 2rjz n iz n — az n — (3zn) 

= -apz n - j3pz n + (a + p)z n + 2/3z n z n H — . (2.30) 

7 + 2 VZn 

Comparison of this system of N ODEs to the system (2.29) yields the following system of N 
algebraic equations (note that - as it were, "miraculously" - the velocities z n have disappeared 
from these equations): 

apz n + Ppzl + [a + p + /« + 2(f) + f^)z n ]f n = M^ + fr£ + /n)A ^ 

7 i £i]z n 

It is now clear that, in order to satisfy this system - identically, i.e. for any values of the 
coordinates z n - one must set either 

case (i) : /„ = (a + bz n )(^ + 2rjz n ) (2.32a) 

implying (see (2.28b)) 

case (i) : / (0) = a 7 , /W = 2ary + 67, / (2) = 26t/, (2.32b) 

or 

case (ii) : /„ = -az„ - /3z 2 + (a + bz n )(j + 2r/z n ) (2.33a) 
implying (see (2.28b)) 

case (ii) : / {0) = 07, = 2m? + 67 - a, / (2) = 26t/ - /3. (2.33b) 

So, in both cases, we hereafter only retain the freedom to assign the 2 constants a, 6 rather 
than the 3 constants /(°), /W, f (2 \ 

Clearly in case (i) we get the following system of N algebraic equations: 

apz n + j3pz^ + [a + p + 2ar/ + 67 + 2(/3 + 26r/)2: n ] [07 + (2arj + 67)2^ + 2br)z 2 ] 

= 2t] [az n + pz 2 n + (a + 6z n ) (7 + 2r/z n )] (a + 6z n ) , (2.34a) 

and likewise in case (ii) 

apz n + (3pz 2 + [p + 2ar] + 67 + 4br]z n ] [07 + (2ar/ + 67 - a)z n + (26r/ - /3)z%] 

= 2r] [-az n - j3zl + (a + 6z n )( 7 + 2r/z„)] (a + 6z n ). (2.34b) 

Hence in case (i) the following set of 4 nonlinear algebraic equations must be satisfied by the 
7 parameters a, b, a, f3, 7, 77, p: 

6r/(/3 + 26r?) = 0, (2.35a) 

(/3 + 26r/)(p + 2ar/ + 267) = 0, (2.35b) 

(a + 2ar/ + 67) (p + 67) + 2a(/3 + 677)7 = 0, (2.35c) 

a(a + p + 67)7 = 0. (2.35d) 

Likewise in case (ii) the following set of 4 nonlinear algebraic equations must be satisfied by the 
7 parameters a, b, a, f3, 7, 7/, p: 

bqifi - 2brj) = 0, (2.36a) 

b[2r](2ar] + 2bj-a + p)- fij] = 0, (2.36b) 

2arjp + 67(4077 + 67 - a + p) = 0, (2.36c) 

a(p + 67)7 = 0. (2.36d) 

Then a trivial if rather tedious computation yields, in case (i) respectively in case (ii), the 
solutions reported in Table 2.1 respectively in Table 2.2 (note that for a = /3 = 77 = these two 
cases coincide, so the corresponding results are only included in Table 2.2). 
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Table 2.1 



# 


a 


b 






7 


V 


p 


1 








* 


* 


* 


* 





2 


* 





* 








* 





3 


* 





—2arj 








* 


* 


4 


* 





—2arj 





* 


* 


2arj 


5 


* 





—2ar] 


* 





* 


—2ar\ 


6 


* 





2ar] 


4ar; 2 /7 


* 


* 


—2ar\ 


7 





* 


—67 





* 





* 


8 


* 


* 


—67 





* 








9 





* 


* 





* 





—67 


10 





* 


—67 


* 


* 





-267 


11 


* 


* 


67 


6 2 7 /a 


* 





-267 


12 





* 


—67 


-2brj 


* 


* 


* 


13 





* 


* 


-2brj 


* 


* 


—67 


14 


* 


* 


* 


-2br] 





* 





15 


* 


* 


—2ar\ 


-2br] 





* 


* 


16 


— &7 


* 


* 


-2br) 


* 


* 


—a — 67 


17 


* 


* 


—2ar\ 


-2brj 


* 


* 


—a — 67 



This table indicates the 17 sets of values to be assigned to the 7 parameters a, b, a, /3, 7, 77, p in order to 
satisfy the 4 equations characterizing case («), see (2.35). Cases with 7 = 77 = are excluded. Asterisks 
indicate that the corresponding parameters can be assigned freely. Note that each of the 7 lines 1, 12, 
13, 14, 15, 16, 17 assigns values to only 3 of the 7 parameters, while each of the other 10 lines assigns 
values to 4 of the 7 parameters. 

Table 2.2 



# 


a 


b 


a 




7 


V 


p 


1 








* 


* 


* 


* 


* 


2 


* 





* 


* 


* 


* 





3 





* 


67 + P 





* 





* 


4 


* 


* 








* 





—67 


5 





* 


bj + p 


2br] 


* 


* 


* 


6 





* 


P 


2br] 





* 


* 


7 


* 


* 


2(27/ 


2br] 


* 


* 


—67 



This table indicates the 7 sets of values to be assigned to the 7 parameters a, b, a, j3, 7, 77, p in order to 
satisfy the 4 equations characterizing case (ii), see (2.36). Cases with 7 = 7; = are excluded. Asterisks 
indicate that the corresponding parameters can be assigned freely. Note that each of the lines 1 and 2 
assigns values to only 2 of the 7 parameters, lines 3, 4 and 6 assign values to 4 of the 7 parameters, and 
lines 5 and 7 assign values to 3 of the 7 parameters. 

We therefore conclude that the many-body models of goldfish type characterized by the 
Newtonian equations of motion (2.30) (or, equivalently, (2.29)) are solvable provided either the 
quantities f n are expressed by (2.32a) and the 7 parameters a, b, a, /3, 7, rj, p are consistent with 
Table 2.1 or the quantities f n are expressed by (2.33a) and the 7 parameters a, 6, a, (3, 7, 77, p 
are consistent with Table 2.2. There are therefore altogether 24 solvable models. Some of these 
models are however not new (in particular, when some parameters vanish); moreover in some 
cases the solution (2.8) with (2.7) of the matrix equation (2.1) is only valid in a limiting sense 
(see Remark 2.2). We leave these issues to whoever will be interested - possibly in specific, 
theoretical or applicative, contexts - in more detailed investigations of anyone of these models; 
our focus in this paper is rather in the unified treatment, and the simultaneous display (see the 
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Table 2.3 








# 


a 


b 


a 




7 


V 


p 


1 





* 


—67 





* 





* 


2 


* 


* 


—67 





* 








3 





* 


* 





* 





—67 


4 





* 


—67 


* 


* 





-267 


5 


* 


* 


67 


6 2 7 /a 


* 





-267 



This table indicates the 5 sets of values to be assigned to the 7 parameters a, b, a, /3, 7, 77, /? in (2.37) 
with (2.38). Asterisks indicate that the corresponding parameters can be assigned freely. Note that in 
every case there are 3 free parameters. 

Table 2.4 

# a b a (3 j 7j p 

1 0*&7 + p0*0 * 

2 * * 0*0 -67 

This table indicates the set of values to be assigned to the 7 parameters a, b, a, ft, 7, 77, p in (2.37) 
with (2.39). Asterisks indicate that the corresponding parameters can be assigned freely: in the second 
line a, b, and 7 are 3 free parameters, in the first line 6, 7 and p are 3 free parameters. 



appendices below), of all of them (except for the elimination, as already mentioned, of the cases 
with 7 = rj = 0). 

It is convenient to write the corresponding equations of motion in two different ways, de- 
pending whether the parameter 77 does or does not vanish. 

If the parameter 77 vanishes, r\ = 0, the equations of motion read as follows: 



v 



-OipZr. 



Ppzl + (a + p)z n + 2f3z n z n + 2(i n + f n ) } 



+ h 



(z n - Zl) 



(2.37) 



with the following assignments corresponding respectively to case {%) and case {it). 
In case {%) 

fn = a7 + b-fz n 



(2.38) 



and the 7 parameters a, b, a, /3, 7, 77, p are restricted according to Table 2.3 (being the relevant 
subcase of Table 2.1). 
In case (ii) 



f n = 07 + (67 -a)z n - /3zl 



(2.39) 



and the 6 parameters a, 6, a, (3, 7, p are restricted according to Table 2.4 (being the relevant 
subcase of Table 2.2). 

If the parameter rj does not vanish, 77 7^ 0, the equations of motion are more conveniently 
written in terms of the dependent variables 



x n (t) = z n (t) + 



27/ '' 



(2.40) 



reading then as follows: 

x 2 



— + A— + j3x n x n + p [x n + A - px n - /3a£) 



N 



{X n + fn) ^2 



(X£ + fe)(Xn + Xj) 
{x n -x£)Xl 



(2.41a) 
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Table 2.5 



# 


a 


b 


a 




7 


V 


p 


1 








* 


* 


* 


* 





2 


* 





* 








* 





3 


* 





—2ar\ 








* 


* 


4 


* 





—2arj 





* 


* 


2ar] 


5 


* 





—2arj 


* 





* 


—2ar\ 


6 


* 





2ar] 


4ar] 2 /j 


* 


* 


—2ar\ 


7 





* 


—67 


-2brj 


* 


* 


* 


8 





* 


* 


-2br] 


* 


* 


—67 


9 


* 


* 


* 


-2brj 





* 





10 


* 


* 


—2arj 


-2brj 





* 


* 


11 


—67 


* 


* 


-2br] 


* 


* 


—a — 67 


12 


* 


* 


—2ar\ 


-2br] 


* 


* 


—a — 67 



This table indicates the 12 sets of values to be assigned to the 7 parameters a, b, a, f3, 7, 77, p in (2.41) 
with (2.42). Asterisks indicate that the corresponding parameters can be assigned freely (with 77 7^ 0) . 
Note that each of the 7 lines 1, 7, 8, 9, 10, 11, 12 assigns values to only 3 of the 7 parameters, while the 
other 5 assign values to 4 of the 7 parameters. 



or equivalently 



( A _ 2/W) — + [(N — 2)/« + p]x n + ( - 2fW + p) ±nXn 

Xn 

+ (N- 2)/(°)/« + [-pp + (N - l)(/«) 2 - 2/(°)/( 2 )]x n 

iV 

+ [-fjp +(N- 2)/«/( 2) K " (/ (2) ) 2 4 + + /„) £ 



(/Co))* 



+ />A 



fc=i 



(i* + / (0) ) , f{2) 



TV 



(X n + / n )Qri + /l 
(z n - 2^) 



+ f^Xk 

(2.41b) 



with 



A 



(2ar] - /?7)7 
4r/ 2 



p = a 



v 



i 2 2 A 

so that a = a — -, 

4 ' 



(2.41c) 



and with the following assignments corresponding respectively to case (i) and case (ii). 
In case (i) 

f n = (2ar] - bj)x n + 2br\x 2 n 



(2.42) 



and the 7 parameters a, b, a, (3, 7, rj, p are restricted according to Table 2.5 (being the relevant 
subcase of Table 2.1). 
In case (ii) 



U = (2aV ^ lh + (2a, -bj-a + ^)x n + (26, - P )xl 



(2.43) 



and the 7 parameters a, b, a, (3, 7, rj, p are restricted according to Table 2.6 (being the relevant 
subcase of Table 2.2): 

These 24 Newtonian equations of motion are exhibited in Appendix A. Their solvable charac- 
ter is of course implied by the fact that the N coordinates z n (t) coincide with the N eigenvalues 
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Table 2.6 



# 


a 


b 


a 




7 


V 


p 


1 








* 


* 


* 


* 


* 


2 


* 





* 


* 


* 


* 





3 





* 


67 + P 


26t? 


* 


* 


* 


4 





* 


P 


2brj 





* 


* 


5 


* 


* 


2an 


2bn 


* 


* 


—67 



This table indicates the 5 sets of values to be assigned to the 7 parameters a, b, a, 0, 7, 77, p in (2.41) 
with (2.43). Asterisks indicate that the corresponding parameters can be assigned freely (with rj 7^ 0). 
Note that each of the first 2 lines assigns values to only 2 of the 7 parameters, while the third and fifth 
lines assign values to 3 of the 7 parameters, and the fourth line assign values to 4 of the 7 parameters. 



of the matrix U(t) (see (2.10a)). To ascertain the behavior of these solutions z n (t) one must 
in each case take account of the restrictions on the parameters characterizing these models, as 
detailed above, which are of course also relevant in order to identify the corresponding evolution 
of the matrix U(t): as implied by inserting in the explicit formula (2.8) with (2.7) - in addition 
to the parameters a, ft, 7, 77, p associated with the solvable iV-body model under consideration - 
the expressions of the initial values U(0) and V(0) = Vq of the matrices U (i) and V(t) in terms 
of the N initial values z n (Q) of the N coordinates and the N initial values z n (0) of the N ve- 
locities. To obtain these expressions it is useful to note that it is possible - and convenient - to 
assume that the diagonalizing matrix R(t) (see (2.10)) is initially just the N x N unit matrix J, 

R(0) = I, (2.44) 

implying initially (see (2.10)) 

U(0) = diag[z n (0)], U nm (0) = eWn(O), (2.45a) 
V = V{0) = Y(0), V nm (0) = 5 nm y n {0) + (1 - 5 nm )Y nm (0). (2.45b) 

The first of these two formulas provides the explicit expression of U(0) in terms of the initial 
data z n (0). 

In the second formula the initial values y n (0) of the diagonal elements of the matrix Vo in 
terms of the initial coordinates z n (0) and velocities z n (Q) of the N particles read 

z n (0) - az n (0) - (3z 2 n (0) 

Vn{0) = t^^GoI (2 ' 45c) 

(see (2.12b)), while the off-diagonal elements Y nm (0) (with n 7^ m) of the matrix Vo are given 
by the ansatz (2.26) yielding 

WO) = {<7n(0)<7 m (0)[i„(0) + /„(0)][i m (0) + / m (0)]} 1/2 , n^m, (2.45d) 
with the quantities g n (0) and / n (0) given by the formulas (see (2.28)) 

^(°) = , nV /» = / (0) + + / (2) 4(0), (2-45e) 

7 + 2rjz n (0) 

with the appropriate assignments of the parameters 7, 77, f(°\ f( 1 ' and characterizing the 
various solvable many-body models, see above (in particular for /^°', f^' and see (2.38) 
or (2.39) or (2.42) or (2.43), as appropriate). 

Note moreover that the dyadic character of the off-diagonal part of the matrix Vo = Y(0), 
see (2.45d), implies a simplification when one must compute functions of this matrix Vq such as 
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those appearing in the explicit expression (2.8) with (2.7) of U(t); this simplification becomes 
particularly significant when the matrix Vq = Y(0) is altogether dyadic, Y nm (0) = v n v m , since 
for any dyadic matrix, say X nm = x n x m , there holds the simple formula 

<p{X) = <p(0)I + ^~^ X, x i = f: xl (2.46) 

k=l 

where tp{x) is any (scalar) function for which the (matrix) expression <p(X) makes good sense. 
This simplification clearly happens iff 

f n = -az n - (3zl, (2.47a) 
as implied by (2.12b) and (2.26) with (2.28), hence in case (i) whenever (see (2.32) and Table 2.1) 

a7 = 0, 2ar] + bj = -a, 2brj = -0, (2.47b) 
and in case (ii) whenever (see (2.33) and Table 2.2) 

a = b = 0. (2.47c) 

Special cases and their (autonomous) isochronous variants. Certain special models among 
those identified above as solvable can be isochronized via the following change of dependent and 
independent variables, 

. , k expfiwt) — 1 ,„ 
z n (t) = exp(kru;t)C„ (r) , r = . ; . (2.48) 

Here the quantities Cn( T ) are assumed to satisfy the Newtonian equations written above, 
see (2.37) with (2.38) or (2.39), of course with the new (complex) independent variable r repla- 
cing the time t; uj is an arbitrary real (for definiteness, positive) constant to which we associate 
the period 

T = —■ (2.49) 

and the number a is adjusted so as to produce, for the dependent variables z n = z n (t) (with the 
real independent variable t interpreted as "time") autonomous equations of motion (the special 
models providing the starting points for the application of this trick being appropriately selected 
to allow such an outcome). Since the application of this trick is by now quite standard (see, for 
instance, Section 2.1 entitled "The trick" of [6]), we dispense here from any detailed discussion 
of this approach and limit ourselves to reporting the results. 

This trick is only applicable to (2.37) with (2.38) in the very special cases with a = p = 
and either 7 = or a = b = (as long as one is only interested in getting autonomous equations 
of motion). Then the assignment a = 1 yields the isochronous equations of motion 



z n = 3iuz n + 2u 2 z n + 2[3z n (z n - iu:z n ) + 2(z n - iuz n ) ) j 



N 



z 7! 



(2.50a) 



Likewise the application of this trick to (2.37) with (2.39), again in the very special cases 
with a = p = and either 7 = 0ora = 6 = 0, and again with the assignment a = 1, yields the 
isochronous equations of motion 

z n = 3iuz n + 2u?z n + 2(3z n (z n - iuz n ) 
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+ 2 (i„ - iujz n - (3zl) 



N 

E 



(2.50b) 



Neither one of these two models is however new: see Examples 4.2.2-6 and 4.2.2-7 in [6]. 

An analogous treatment applied, mutatis mutandis, to (2.41) with (2.42) in the special case 
with p = 0, 2arj = 67 and 2ar\ = fij (implying A = 0), yields (again via the assignment a = 1) 
the isochronous equations of motion 



\ujx n +lo x n + — + f3x n (x r , 



iuix n ) + (x n - \ojx n + 2bnx n ) 



N 

E 



(x£ — iux£ + 2brjxf){x n + xi) 
(x n - xt)xi 



(2.50c) 



Likewise an analogous treatment applied to (2.41) with (2.43) in the special case with p = 
and either 2ar/ = (3j (implying A = 0) and a = —2ar\ + 67 or 7 = (implying A = 0) and 
a = 2ar], yields (again via the assignment a = 1) the same isochronous equations of motion (up 
to the, merely notational, replacement of 2br] with 2brj — /3). 

While finally this treatment applied, with a = —1, to (2.41) with (2.43) in the special case 
with b = (3 = p = and a = 2ar\ yields the isochronous equations of motion 



lOJXr. 



2 yj 3*71 • 

u x n -\ h 07 h 0710; 



N 



(x n + lUJX n + 07) ^ 



(xi + iujxe + aj)(x n + xj) 

(x n - Xi)X£ 



(2.50d) 



2.3 A related class of solvable many-body models 

In this subsection we consider the Newtonian equations of motion that obtain by identifying 
the dependent variables of the models discussed above as the N zeros of a monic (time- 
dependent) polynomial of degree N, and by then focussing on the time-evolution of the N 
coefficients of this polynomial. It is again convenient to treat separately the two cases with 
r] = and with 77 7^ 0. 

In the rj = case the starting point are the equations of motion (2.37) with (2.38) or (2.39). 
We then introduce the time-dependent (monic) polynomial ip(z, t) whose zeros are the N eigen- 
values z n {t) of the N x N matrix U(t): 



dei[zl 

N 



A' 



,N-ml 



(2.51a) 
(2.51b) 



n=l 



m=l 



The last of these formulas introduces the N coefficients c m = c m (t) of the monic polynomial 
tp(z,t); of course it implies that these coefficients are related to the zeros z n (t) as follows: 



N 



Cl 



C-2 



n=l 



N 

E 

n,m=l]n>m 



Z n Z r . 



(2.51c) 



and so on. 

The fact that the initial- value problem associated with the time evolution of the coordi- 
nates z n can be solved by algebraic operations implies that the same solvable character can be 
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attributed to the time evolution of the monic polynomial vp(z, t) and of the N coefficients c m (t). 
The procedure to obtain the equations of motion satisfied by the N coefficients c m (t) from 
the ./V equations of motion satisfied by the N zeros is tedious but standard; a key role in this 
development are the identities reported, for instance, in Appendix A of [6] (but note that there 
are two misprints in these formulas: in equation (A. 8k) the term (N + 1) inside the square 
brackets should instead read (N — 3); in equation (A. 81) the term N 2 inside the square brackets 
should instead read N(N — 2) - these misprints have been corrected in the recent paperback 
version of this monograph [6]). Here we limit our presentation to reporting the final result. 

The equation characterizing the time evolution of the monic polynomial ip(z,t) implied by 
the Newtonian equations of motion (2.37) with (2.28b) reads as follows: 

At ~ 2(/ (0) + f W z + + (pV>+p<»z)ilH 

+ (4n 4 x w gf, 2 + 4 3 w ^ 4 )^ 

+ + q{ 1] z + qfz 2 + gfVty, + + + z 2 )^ = 0, (2.52) 



with 



(°) =- a - p + 2(N- 2/ 2 ) ci, p (1) =2[-/3 + (7V-2)/( 2 )]; (2.53a) 



P 

4 0) = (/ (0) ) 2 > q? = 2/ (0) / (1) , q? ] = 2/ (0) / (2) + (/ (1) ) 2 , 



4 3) =2 / «/( 2 ), 4 4) = (/ (2) ) 2 ; (2.53b) 
,(0) = _ 2(A r_ 1)/ (0) / (i) + 2/ (0) / (2) Ci; 

q 0) = - ap - 2 {N - 2)/(°)/( 2 ) + 2/( 1 )/( 2 ) Cl - 2{N - 1) (fW) 2 , 
^ 2) = -/3p-2(2A-3)/( 1 )/( 2 ) + 2( / ( 2 )) 2 c 1 , 

q[ 3) = -2(N-2){fW) 2 ; (2.53c) 

= Nap - 27V/(°)/( 2 ) + N(N - l)(/«) 2 

- \p P + 2(AT - l)/«/( 2 )] Cl + 2(/3 + f^) Cl + 2(/( 2 )) 2 c 2 , 

= NPp + 2A(iV - 2)/( 1 )/( 2 ) - 2(N - 1)(/< 2 >) 2 ci, 
q {2) =N(N-3){f^) 2 , (2.53d) 

where of course the quantities /^°^ , should be expressed in terms of the other parameters 

as implied by (2.28b) with (2.38) or (2.39), as the case may be. Note that, while this equa- 
tion, (2.52), satisfied by the function ip(z, t) (where of course subscripted variables denote partial 
differentiations) might seem a linear PDE, it is in fact a nonlinear functional equation, because 
some of its coefficients, see (2.53), depend on the quantities c\ and c 2 which themselves depend 
on ip, indeed clearly (see (2.51)) 

cm = c m (t)= ^ {N _l y (2.54a) 

where we used the shorthand notation t/j^\z, t) to denote the j-th partial derivative with respect 
to the variable z of ip(z,t), 

i& Mm **M, , = 1,2,.... (2.54b) 

Likewise, the equation characterizing the time evolution of the monic polynomial 4> { x i t) 
implied by the Newtonian equations of motion (2.41) via the following assignment (analogous 
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to (2.51b)), 

N N 

4>{x, t) = l\[x- x n (t)\ =x N + J2 [c m (t)x N - m ] , (2.55) 

n=l m=l 

reads 

+ q W x + g <V + g (3) x 3 + ^4)^ (2 _ 56) 

+ f + ^ + q ( o ] x + sf^l = 0, (2.57) 



X 



now with 



= 2/ (0) _ A) (0) = _ p + Nf (D _ /(2) Cl + jW^V-l _ PAT 

CiV Cat 



P 



(1) 



92 



(0) 



2 (iV_i) / (2)_ /3 . (2.5 8a ) 
4 3) =2 / «/( 2 ), 4 4) = (/ (2) ) 2 ; (2.58b) 

gC-D = .(/(O)) 2 , ,(0) = pA _ ^(0)^1) + /( 0) /(2 ) Ci + /(0 ) PAT _ (/ (0) ) 2^ ) 

CAT Cat 

g« = -p M - 2(iV - l)/(°)/( 2 ) + /«/( 2 ) Cl - (iV - l)(/«) 2 + - /W/W^i, 

Cat Cat 



(2) _ _ {3N _ 4)/ (D / (2) + (/ ( 2)) 2 ci + /(2) CAT _ /(0 ) fi2 )?N-l 



1i 

CAT Cjy 

g (3) = _(2Ar_ 3 )(/( 2 )) 2 ; (2.58c) 
^ = (A-2/W)^ + (/(0)) 2 ^, 

CAT CAT 

g?> = -/3/> Cl + ^ + (/ (2) C1 - AT/(D)^ + f (°) {Nf (D _ /(2) Cl) ^ ; 

CA" ' c n 

« = iV/3p + iV(iV - 2)/W /( 2 ) - (iV - l)(/ (2) ) 2 Cl - iV/( 2 )^ + /(°) /( 2 )^ 



96 

CAT CAT 

q W =N(N-2){fM) 2 . (2.58d) 

The equations of motion of Newtonian type satisfied by the N coefficients c m (t) which obtain 
from (2.52) hence correspond to the Newtonian equations of motion (2.37) read as follows: 

cm - 2{N - m + l)/ (0) c m _i + [ - 2{N - m)/W + p {0) ] c m 

+ [- 2 (N -m- l)/^ + p (i)] c m+1 + (AT - m + 2)(N - m + l)4 0) c m _ 2 
+ (N - m + 1) [(N - m)q { 2 1) + q^] c m _i 
+ {(N - m)[(N - m - l)gf + q?) + q^}c m 
+ { (AT - m - 1) [(N - m - 2)qf ] + g{ 2) ] + }c m+1 
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+ { (AT - m - 2) [(A - m - 3)q ( 2 4) + q^] + q^ }c m+2 = 0, (2.59) 

where c n vanishes for n < and for n > A while Co = 1 (see (2.51)), and of course the 
coefficients p^ and qip are defined by (2.53). Again, this system of ODEs might seem linear, 
but it is in fact nonlinear because some of its coefficients depend on the dependent variables c\ 
and C2, see (2.53). The more explicit version of these equations of motion that obtain by 
expressing the various coefficients in terms of the free parameters are listed in Appendix B. 
They are of course just as solvable as the Newtonian equations of motion satisfied by the A 
coordinates x n (t), see Appendix A, to which they correspond via (2.51); and in particular 
whenever the parameter p is imaginary and the parameter a is real and negative they are 
asymptotically isochronous with period T (see Remark 2.3). 

Likewise, the equations of motion of Newtonian type satisfied by the A coefficients c m (t) 
which obtain from (2.57) hence correspond to the Newtonian equations of motion (2.41) read 
as follows: 

c m + [p { ' 1] -2{N-m + c m -i + [-2{N - m)/« + p^] c m 

+ [-2(N -m- 1)/^ + p W] dm+1 + (A - m + 2) [(N - m + l)q!® + q^] c m -2 

+ { ( N - m + 1) [( A - m)4 1} + q?] + </Jf 1} }c m -i 

+ {(A - m)[(N - m - l)q® + q®] + q^}c m 

+ { (AT - m - 1) [(A - m - 2)qf ] + qf ] ] + q^ ) }c m+1 

+ { (AT - m - 2) [( A - m - 3)q^ + q[ 3) ] + q^ }c m+2 = 0, (2.60) 

where of course again c n vanishes for n < and for n > N while cq = 1 (see (2.55)) and of 
course the coefficients p^ and qi are now defined by (2.58). Again, this system of ODEs 
might seem linear, but it is in fact nonlinear because some of its coefficients depend on the 
dependent variables ci, cn-i and cat, see (2.58). The more explicit version of these equations 
of motion that obtains by expressing the various coefficients in terms of the free parameters are 
listed in Appendix B. They are of course just as solvable as the Newtonian equations of motion 
satisfied by the A coordinates x n (t), see Appendix A, to which they correspond via (2.55); and 
in particular whenever the parameter p is imaginary and the parameter a is real and negative 
they are asymptotically isochronous with period T (see Remark 2.3). 

Special cases and their (autonomous) isochronous variants. Certain special models among 
those identified above (in this subsection) as solvable can be isochronized by an analogous trick 
to that employed at the end of the preceding Subsection 2.2. One route to this end takes as 
starting point the isochronized systems of Newtonian equations of motion (2.50) and applies to 
them the same procedure employed above to obtain the equations of motion (2.59) with (2.53) 
and (2.60) with (2.58). An equivalent procedure is to apply to certain special subcases of these 
systems of ODEs, (2.59) with (2.53) and (2.60) with (2.58), the following change of dependent 
and independent variables: 

exp (iuit) — 1 

Cm{t) = exp(i amojt) \m (t) ,t = ; . (2.61) 

Here the quantities \ n (r) are assumed to satisfy the systems of ODEs written above, see (2.59) 
with (2.53) and (2.60) with (2.58), of course with the new (complex) independent variable r 
replacing the time t; uj is an arbitrary real (for definiteness, positive) constant to which we 
associate the period T, see (2.49); and the number a is adjusted so as to produce autonomous 
ODEs for the new dependent variables c m = c m (t) (with the real independent variable t inter- 
preted as "time": the special models providing the starting points for the application of this 
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trick being appropriately selected in order to allow such an outcome). Since the application of 
this trick is quite standard, we dispense here from any detailed discussion of this approach and 
limit ourselves to reporting the results. 

The ODEs that follow from (2.59) (with a = p = rj = 0, /(°) = /M = = and a = 1) 
read as follows: 



dm = i(2m + \)ujc m — [2/3(ci — iwci) — m(m + l)w 2 ]c r , 
+ 2/3c m +i - i2/3(m + l)wc m+ i; 



(2.62a) 



those that follow from (2.59) (with a = p = 77 = 0, /(°) = /W = 0, = -/3 and a = 1) read 
instead as follows 



[i(2m + l)w — 2/3ci]c m + [m(m + l)u> 2 + i2m/3ojci — 20 2 C2\c m 

+ 2mc m+ i + 2m[/3 2 ci - i(m + l)cj]c m+ i — (m— l)(m + 2)/3 2 c m+2 . 



(2.62b) 



Neither one of these two isochronous many-body problems is new. 

The analogous results that follows from (2.60) are instead generally new. There are then two 
sets of cases. The first set of isochronous models obtain from the assignment a = 1 and read 



-i(N - 2m - l)u + / (2) ci + 



cn 



+ {P- 2mfW) 



Cm+l 



+ 



+ 



-m(N - m - l)w 2 + (N - m)w/ (2) ci - —(imu + / (2) ci + 0(iuQ - ci)) 

cn 



-UN - 2m) (m + l)w/ (2) - Um + + m(f (2) ) 2 Cl + (m + l)/ (2) — 

cn 



Cm+l 

(2.63a) 



-m(m + 2)(/( 2 )) 2 c m+2 
with the following restriction on the parameters: 

p = A = /(°) = /W = 0. (2.63b) 
The second set of isochronous models obtain from the assignment a = — 1 and read 

c m = (2N - 2m + l)Ac m _i + 



i(iV-2m + l)u; + ^-A^ 

Cat C7V 



(N - m)(N - m + 2)A 2 c m _ 2 + 



Um - 1)(2N - 2m + l)w 



+ IN - m + 1)(JV - m - 1)A ( A^=^ 

V CAT 



iiV^ - k JL] 

CN J 



Cm-1 



+ 



-m{JS — m + ljw + lmw I A 

cat Cat 



with the following restrictions on the parameters: 

b = /3 = p = /« = /( 2 ) =0, a = 2ar?, /<°> = A 



07. 



(2.64a) 



(2.64b) 



3 Outlook 

Results analogous, but somewhat more general, than those reported in this paper can be obtained 
by an analogous treatment based on a somewhat more general - but still solvable - system of 
two N x N matrix ODEs than (2.1), such as, for instance, 



U = aU + (3U 2 + jV + V (UV + VU), V = p + pV + p 2 V 2 , 



(3.1) 
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which contains the 2 additional scalar constants po and p2 (and clearly reduces to (2.1) for 
Po = p2 = 0). These developments will be reported in subsequent papers. 

Finally, let us recall that Diophantine findings can be obtained from a nonlinear autonomous 
isochronous dynamical system by investigating its behavior in the infinitesimal vicinity of its 
equilibria. The relevant equations of motion become then generally linear, but they of course 
retain the properties to be autonomous and isochronous. For a system of linear autonomous 
ODEs, the property of isochrony implies that all the eigenvalues of the matrix of its coefficients 
are integer numbers (up to a common rescaling factor). When the linear system describes 
the behavior of a nonlinear autonomous system in the infinitesimal vicinity of its equilibria, 
these matrices can generally be explicitly computed in terms of the values at equilibrium of the 
dependent variables of the original, nonlinear model. In this manner nontrivial Diophantine 
findings and conjectures have been discovered and proposed: see for instance the review of such 
developments in Appendix C (entitled "Diophantine findings and conjectures") of [6]. Analogous 
results obtained by applying this approach to the isochronous systems of autonomous nonlinear 
ODEs introduced above - and in subsequent papers - will be reported if they turn out to be 
novel and interesting. 



A First appendix 

In this appendix we list the 24 Newtonian equations of motion whose solvable character has 
been demonstrated in this paper. In each case the parameters they feature (such as a, b, a, 
P, 7, r\, p, as the case may be) are arbitrary constants; the (assigned) values of the other 
ones of these parameters (which also characterize the time-evolution of the solutions of these 
equations, see (2.8)), are also reported. Let us emphasize that if the parameter p is an imaginary 
number and the parameter a is real and negative, the corresponding many-body problem is 
asymptotically isochronous with period T, see Remark 2.3; and that isochronous many-body 
models are characterized by the 4 Newtonian equations of motion (2.50) displayed at the end of 
Subsection 2.2. Let us also mention again that the equations of motion reported below are not 
all new; in particular not new are clearly those whose corresponding equations of motion in the 
following Appendix B are linear. 

rj = 0, case (i) (5 models, corresponding to Table 2.3): 

(1) a = P = 0, a = -67, /„ = -az n : 

n fz-az\ 

z n = -apz n + (a + p)z n + 2(z n - az n ) ^ I - - - j ; (A. la) 

(2) a = -67, P = p = 0, f n = -(a/b)a - az n : 



, / . act \ f Z£ — aa/b — aze \ 
z n = -az n + 2 [z n - — - az n J ^ I — — j; (A.lb) 

(3) a = ft = 0, p = -67, /„ = -pz n : 



N 



z n = -apz n + (a + p)z n + 2(z n - pz n ) ^ f _ j ; (A.lc) 

e=i, i^n ^ Zn Zi ' 

(4) a = 0, a = -67, p = -2bj = 2a, f n = -az n : 

N ( z -az\ 

z n = -2a 2 z n -2af3zl + 3az n + 2Pz n z n + 2(z n - az n ) V" ( ); (A.ld) 
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(5) a = 67, (3 = b 2 j/a = ba/a, p = —267 = —2a, f n = aa/b + az n : 



2b 



2b 



z n = 2a 2 z n H c?z\ + — az n z n - az n 



N 

+ 2 [z n + + az„j ^2 



Z£ + aa/b + az£ 
z n - zt 



7] = 0, case (m) (2 models, corresponding to Table 2.4): 
(1) a = ft = 0, a = 67 + p, /„ = -pz n : 



N 



-(67 + p)pz n + (67 + 2p)i n + 2(z n - pz n ) ^2 



--l,t±n 



zt - pz t 

Z n ~ Zt 



(2) a = P = 0, p = -67, / n = 07 + bjz n : 



N , 

-b^z n + 2(z n + 07 + &7z n ) ^ ( 

t=l,t^n ^ 



il + 07 + 67^ \ 

Zn~ Zt J 



7] 7^ 0, case (i) (12 models, corresponding to Table 2.5): 
(1) a = b = p = 0, /„ = 0, A = (2ar, - ^7)7/ (V), M = <* - 



JV 



X n — -|- A -|- pX n %n X n / 

'i ' 1 • & / 

i=l, tytn 



x e (x n + Xt) 



(x n - xt)xt\ ' 
(2) b = = 7 = p = 0, /„ = 2ar?x„, A = 0, p = a: 



JV 



x n = — + (x n + 2ar]x n ) ^2 

Xn t=l,l+n 



(±e + 2ar\xt)(x n + x e ) 



(x n -xt)xt 

(3) b = P = 7 = 0, a = —2ar], f n = —ax n , A = 0, p = a: 



A? 



(x£ - ax e )(x n + 



— "I - p{x n QXn) ~\~ (x n ttX^) ^ ^ 

x?1 e=i,e^n 



(4) 6 = /3 = 0, a = -2ar?, p = -a, / n = -ax n , A = 07/ (2?]), p = a: 



■■ _< ajxn _ . ,<^Y_ 
x„ 2i] x n \ 2r\ 



N 



+ (x n ~ ax n ) ^ 

1=1, tyn 



(xe - axt)(x n + Xt) 



(x n ~ Xt)Xt 

(5) b = 7 = 0, a = p = -2ar], f n = -ax n , A = 0, p = a: 
-\- Px n x n ~\~ ol (x n ax n /3x n ) 



X I, 



N 



+ (x n OLXn) ^ ^ 

1=1, t^n 



(xt - axt)(x n + xe) 
(x n - xt)xt 



(A.le) 



(A.2a) 



(A.2b) 



(A.3a) 



(A.3b) 



(A.3c) 



(A.3d) 



(A.3e) 
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(6) b = 0, a = 2ar], f3 = Aarj 2 /^ = 2rja/^, p = —a, f n = az n , A = 0, p = —a: 



x\ 2arj . 



x n — -\- X n X n OL I x n -\- OiX n x n 



%n 7 



2arj 



N 



+ (x n + ax n ) 

1=1, tyn 



7 



(±e + axi)(x n + xt) 



(x n - x e )x e 

(7) a = 0, a = -bj, (3 = -26??, f n = ax n - (3x n , A = 0, p = -a: 



(A.3f) 



x n — -\- {3x n x n -\- p (x n -|- ax n -\- /3x n ) 

Xn 

r (±e + axe - fixf) (x n + x£) 



+ (x n + ax n - I3x 2 n ) ^ 

1=1, tyn 



(x n - Xi)x£ 



(A.3g) 



(8) a = 0, = -2bq, p = -67, /n = -&7^n + ^7]x 2 n , A = (a + 67)7/(2??), M = « + 267: 



x ra — + 



x 2 n (a + b"/)jx n 



- 2bi]x n x n - 67 



Xn 

27] X n 

+ (x„ - 67x n + 2br]x n ) 2^ 

e=i,e^n 



x n + ^ 67 ^ 7 - (a + 26 7 )x„ + 2br)x 2 n 

2,17 



(x n - x e )x e 



(A.3h) 



(9) 7 = p = 0, (3 = 2b7], f n = 2a7jx n + 2brjx 2 , A = 0, p = a: 
x 2 

x n = — - 2b7]X n X n + (i n + 2a7]X n + 2b7]X n ) 



N 

* E 

1=1, t^n 



(±£ + 2a7]X£ + 2&7?xf ) (x n + Xi) 
(x n - Xl)Xl 



(A.3i) 



(10) 7 = 0, a = — 2a?7, /3 = — 26?y, / n = — ax n — A = 0, ji = a: 



+ p (x n - ax n - f3x 2 n ) 

(xe - ax e - Pxj) (x n + xe) 



N 



+ (x n ~ ax n ~ Px 2 n ) 

1=1, tyn 



(x n - x e )x£ 



(A.3j) 



(11) a = -67, (3 = -2b7], p= -a- 67, f n = -(2r? + l)6 7 x n + 2bi 1 x 2 n , A = (a + 67)7/(2??), 
p = a + 267: 



_ in (a + 67)7 x n 

x n — ~r 

Xn ^] Xn 



2b7]x n x n - (a + 67) 
+ [x n - (2r? + l)6 7 x n + 26r?x^] 



x n + ^ a ± &7 ^ 7 - (a + 26 7 )x n + 2b7]x 2 n 



2r? 



AT 

■< E 

£=l,t^n 



[xe - (2t] + l)bjxe + 26r?x|] (x n + x t ) ^ _ 

(^n ■ 



(A.3k) 



(12) a = -2ar?, /? = -26r?, p = -a - 67, / n = (2ar? - 6 7 )x„ + 26r?x2 , A = (-2ar? + 67)7/(21?), 
p = — 2ar? + 267: 



.. _ xj | (67 - 2ar?)7 x w nu . 

X n — -f- ZOT}X n X n 
X n Z7] X n 
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- (a + 67) x n + 2a7 l)~f + (2ar? _ 2b 1 ) Xn + 2br]x 2 n 

+ [x n + (2arj - 67)x n + 26ryx^] 

f [xg + (2an - bj)xg + 2brjxf\ (x n + x £ ) 1 



* E 



(x n - x^ 



(A.31) 



77 7^ 0, case (m) (5 models, corresponding to Table 2.6): 

(1) a = b = 0, /„ = A - fxx n - 0x1, A = (2arj - ^7)7 / '(4t? 2 ) , p = a- ^7/77: 

in = — + A— + f3x n x n + p (x n + A - /ix n - /Jx 2 ) 

■En -^n 



TV 



+ (x n + a - - ^ic^) 



(±£ + A - pxg - f3xj) (x n + xg) 



(x n - xg)xg 

(2) b = p = 0, f n = A + (2ar? - /j)x„ - /3x 2 , A = (2ar/ - /3 7 )7/(4r ? 2 ), p = a- /3j/rj: 

+ A^ + /3x n x n + [x n + A + (2ar/ - /i)x n - /3X 2 ] 
[x^ + A + (2a?7 - ^)x^ - /3xj] (x n + xg) 1 

(X„ - X£)X£ J ' 

(3) a = 0, a = 67 + /?, (3 = 2brj, f n = A - /ox n , A = (p - 67)7/(277), p = p - bj: 



(A.4a) 



X?2 X77, 

TV 

E 

g=l,i^n 



(A.4b) 



— + A— + 2br\x n x n + p [x n + A - (p - &7)x n - 26?7x^ 
x n x ra 



TV 



+ (x„ + a - px n ) 

g=i,e^ n 



(xg + A - pxg)(x n + x^) 



(x n - x^)x^ 

(4) a = 7 = 0, a = p, j3 = 2br], f n = -ax n , A = 0, p = a: 
■A 

— +/3x 



(A.4c) 



- 1 • 



TV 



+ (x n - ax n ) y j 

g=l t g^n 



(xg - axg)(x n + xg) 



(x n - Xg)Xg 

(5) a = 2ar/, P = 2 br], p = -67, f n = A + &7x n , A = (2a?7 - ,57)7/(27?), p = 2arj- 267: 

— + A— + 2br\x n x n - 67 [x n + A - (2a7/ - (3j)x n - 2br]x^\ 
x n x n 



(A.4d) 



TV 



+ (x n + a + 67x n ) 

g=l/^n 



(xg + A + bjxg)(x n + x^) 
(x n - x^)x^ 



(A.4e) 



B Second appendix 

In this appendix we list the second series of 24 Newtonian equations of motion whose solvable 
character has been demonstrated in this paper; they correspond to those reported in Appendix A 
via the transformation among the tV zeros z n and the tV coefficients c m of a monic polynomial, 
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see (2.51) and (2.55). In each case the parameters they feature (such as a, b, a, /3, 7, 77, p, as the 
case may be) are arbitrary constants; the (assigned) values of the other ones of these parameters 
(which also characterize the time-evolution of the solutions of these equations, see (2.8)), are also 
reported. Let us emphasize that if the parameter p is an imaginary number and the parameter a 
is real and negative, the corresponding many-body problem is asymptotically isochronous with 
period T, see Remark 2.3; and that isochronous many-body models are characterized by the 4 
Newtonian equations of motion (2.62), (2.63) and (2.64) displayed at the end of Subsection 2.3. 
Let us also mention again that the equations of motion reported below are not all new; in 
particular all those that are linear are of course not new. 

Let us recall that it is always assumed that c n = for n < and for n > N, and cq = 1. 

77 = 0, case (i) (5 models, corresponding to Table 2.3): 

(1) a = p = 0, a= -67: 

Cm + [(1 - 2m)a - p]c m + ma[(m - l)a + p]c m = 0; (B.la) 

(2) p = p = 0, a = -b T . 

c m - 2(N -m- l)a7c m _i + (2m - l)&7c m + (N - m + 2)(N - m + l)a 2 j 2 c m - 2 

- 2(N - m + l)(m - l)ab-f 2 c m -i + m(m - 1)&V = 0; (B.lb) 

(3) a = p = 0, p = -67: 

c m + [(1 - 2m)p - a]c m + mp[(m - l)p + a]c m = 0; (B.lc) 

(4) a = 0, a = —67, p = 2a: 

c m - (2m + l)ac m - 2/3c m +i + [m(m + l)a 2 - 2a(3ci + 2/3ci]c m 

+ 2(m + l)a(3c m +i = 0; (B.ld) 

(5) a = 67, P = b 2 j/a, p = -2bj: 

c m - 2(N -m + l)ajc m -i + (2m - l)67c m - ^-^-c m +i 

a 

+ (N - m + 2)(N - m + l)a 2 -f 2 c m -2 - 2{m -l)(N-m+ l)abj 2 c m -i 



+ 



mym — 0)0 7 H ci H c\ 



2 
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(m + 1) -c m +i = 0. (B.lc) 



77 = 0, case (ii) (2 models, corresponding to Table 2.4): 

(1) a = p = 0, p = -67: 

Cm ~ 2(N - m + l)ajc m -i + (2m - l)&7C m + (N - m + 2)(N - m + l)a 2 -y 2 c m - 2 

- 2(m -l)(N-m+ l)abj 2 c m -i + m(m - l)6 2 7 2 c m = 0; (B.2a) 

(2) a = P = 0, a = 67 + p : 

c m - (2mp + 6 7 )c m + [mb-fp - (N 2 - m 2 ) p 2 ] c m = 0. (B.2b) 

7/ 7^ 0, case (i) (12 models, corresponding to Table 2.5): 
(1) a = 6 = p = 0, A = l(2a V - /3 7 h]/(2r/) 2 , M = a - pj/ V : 

c m ~ Ac m _i - —c m - (3c m+ i + \—c m -i + {5c\c m = 0; (B.3a) 

CN C N 
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(2) b = p = 7 = p = 0, A = 0, p = a: 

c N 



(2N - Am)ar] 



cn 



c m - 2marj ( 2ar\ + — ) c m = 0; 

cn J 



(3) 6 = ft = 7 = 0, a = —2am A = 0, p = a: 

cn 



Cm 



(N-2m)a- p 



cn 



Cm 



m{N — m)a + (N — m)ap — ma 



cn 
cn 



(B.3b) 



0; (B.3c) 



(4) b = ft = 0, a = —2ar], p = 2arj = —a, A = —07, p = a: 

CN 



Cm ~t~ 07 c m— 1 



+ 



(AT - 2m + l)a + 

Cat 

,2 , r>/A r ™A„2 



-(iV - m)a z + 2(N - m)a z + ma 



c m + a<y 
cn 



CN 



(N -m - l)a 



0; 



CN 
CN 



Cm— 1 



(B.3d) 



(5) b = 7 = 0, a = p = —2am A = 0, p = a: 

cn 



Cm ~\~ 



(N -2m- l)a 



cn 



cn 



(m + l){N — m)a + ma \- a(3c\ — (3c\ 

cn 



+ {m + l)a/3c m +i = 0; 
(6) b = 0, a = 2am f3 = Aarj 2 /^, p = —2am A = 0, p = —2am 

CN 



(B.3e) 



Cm 



2(N - 2m + l)ar] + 



CN 



— 2ar\ 



Aarj 2 . 

Cm C m -\-l 

1 

cn ^arj 2 



2/i 

2(m — l)(iV — 777)077 + 777— — — ci ci 

CAT 7 7 



(B.3f) 

g(j2„3 

c m - (ttt, + 1) c m +i = 0; 

7 



(7) a = 0, a = -67, Q = -2bm A = 0, p 

cn 



-a: 



Cm. 



(N - 2m)a + p + /3ci + 



(2777 + l)/3c m+ l 



CAT 



m(iV — m)a — (iV — m)ap + (iV — m)aftci + (/3ci + ma) h /3pci — /3ci 

Cat 



Cm 



- <^ (m + l)/3/9 - [2Ar _ (3 m + 5)tv + 4m + 4]a/3 + m/Tci + (m + l)/3— } c m +i 

cat 



- m(m + 2)/3 2 c m +2 = 0; 
(8) a = 0, = -2677, p = -67, A = (a + 67)7/(27/), /i = a + 26 7 : 



(B.3g) 



ACm-l + 



(N - 2m + 1)67 - 267/ci 



cn 
cn 



C m + 4(777 - l)c m+ i 



+A 



■(N-m + 1)67 + 



CAT 
CJV 



Cm— 1 



— (N — m)ba^ + 2br)C\ + (77767 — 26t/ci 



(m + 4)(JV - m)&V + 2(2jV - 2m + 1)6^77701 
c/v 



cjv 



+2 



(m + l)(iY - 2m + 1)6^77/ - 2mb 2 ifci - (m + 1)677 



cjv 
/ — 

CJV 



Cm+1 



+ 4m(m + 2)c m+2 = 0; 



(B.3h) 
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(9) 7 = p = 0, /3 = -2brj, A = 0, p = a: 

c N 



Cm. 



2(N - 2m) + 2br]d + 



CAT 



c m + 2(2m + l)br]c m+ i 



+ 



—Am(N — m)a 2 rj 2 + 4(N — m)abrj z c\ + /3ci + 2{brjC\ — mar]) 



cn 



- 2 



2(m + 1)(N - 2m)abrj z + 2mb 2 rfc 1 + (m + 1)677 



Cat 

/ 

Cat 



Cm+l 



+ 4m(m + 2)c m+2 = 0; 
(10) 7 = 0, a = -2ar/, /3 = -26??, A = 0, /j = a: 

CAT 



(B.3i) 



Cm ~l~ 



(iV - 2m)a - p + ^ci - 



CAT 



c m - (2m + l)/3c m +i 



m(N — 7x1)0? + (iV — m)ap — (N — m)afic\ + fipci — f3d\ + (/3ci — ma) 



cn 
cn 



(m + 1)(N- 2m)af3 - (m + l)f3p + m(3 2 a - (m + l)/3 



CjV 

cn 



Cm+l 



+ m(m + 2)c m+ 2 = 0; 
(11) a = —67, /3 = — 2br], p = —a — 67, A = (a + 67)7/(2??), 77 = a + 267: 



(B.3j) 



Cm - Acm-i + < [2(N - 2m)r) + N - 2m + l]bj + a - 2br]d - 



cn 



+ 2(2m + l)c m+ i + A 



cat 



l2„,2/o^ 1 i\2 



- (N -m + l)(a + &7) 



Cm— 1 



m(N - m)b 2 Y{2v + 1) - (N ~ m)(a + bj)(a + 267) + 2(N - m)b 2 jrj(2r] + l)ci 



+ 2brj(a + 67)01 + 2677C1 — [2brjC\ + 77167(277 + 1" 



■ CAT 
CAf 



+ 2 



(m + l)(iV - 2m)b 2 jr](2r] + 1) + (m + 1)677(0; + 67) - 2mb 2 r] 2 c 1 



— (m + 1)6?7 



CjV 

cn 



Cm+i + m(m + 2)c m +2 = 0; 



(B.3k) 



(12) a = — 2ar/, /3 = —2bm p = —a — 67 = 2a?? — 67, A = (67 — 2077)7/ (2r?), p = 267 — 2077: 



c — Ac m -i — 



(N - 2m + l)(2ar/ - 67) + 26r7Ci + 



cat 



+ 2(2m + l)677Cm+i + A 



(N -m- l)(2a?7 - 67) + 



CN 



Cm— I 



+ <^ - m(N - m)(2ar/ - bj) 2 + 2(N -m + l)(2ar/ - 67)67701 



+ (N — 771) (2a7i - 67)(2ar? - 267) - 2677C1 + [267701 - 771(2017 - 67)]— f>c m 



(m + l)(iV - 2m + l)(2a7? - 67)677 + 2mb 2 r] 2 c 1 + (m + 1)677 



+ 4m(m + 2)&Vc m+2 = 0. 



cn 
1 — 

CAT 



Cm+l 
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77 7^ 0, case (ii) (5 models, corresponding to Table 2.6): 
(1) a = b = 0, A = [(2a V - ^) 1 ]/(2i 1 ) 2 , p = a- fr/r,: 



c m - (2TV - 2m + l)Ac m -i + 



(TV - 2m)p - p + (3d + A 



Cat Cat 



- (2m + l)/3c m+ i + (iV - m + 2)(TV - m)X 2 c m - 2 - 



(TV -2m){N -m + l)Xp 



(TV - m + l)pX + (N-m + l)/3Aci + A— - X 2 ^± 

c N c N 



Cm— I 



+ 



m(N - m)(2/3A - p 2 ) + (/3ci - mp) (A 



CJV-1 Cat 



V Cat Cat 



-(TV - m)pp + (TV — m)/3pc 1 - Ppc 1 + /3ci 



(m + 1)(TV - 2m)(3p + (m + + m/3 2 ci - (m + l)/3 



CN 
CN 



(N-m- 2)/3A 



CJV-1 
CN 



c m +i + m(m + 2)0 2 c m+2 = 0; 



(2) 6 = p = 0, A = [(2m? - ^) 7 }/(2rj) 2 , p = a — fa/rj: 



Cm - (2TV - 2m + l)Ac m -i + 



(/V - 2m) (p - 2ar]) + 0a + A 



CAT-l Cat 



CAT 



(2m + l)/3c m+ i + (iV - m)(TV - m + 2)A 2 c m -2 



+ 



(TV - 2m) (TV - m + l)X(2ar) — p,) — (N — m+ l)/3Aci 

c m _i + i m(TV - m) [2/3 A - (2ar/ - p) 2 ] 



+ (A r_ m ) A ^-(TV-m)A 2 ^ 

CAT Cat 



- (TV - m)(2arj - p)f3ci + /3ci + [fid + m(2ar? - pi)] A 



CJV-1 _ c^v 
CN C N 

(m + 1)(TV - 2m)/3// + (m + l)/3p + m,/3 2 ci 

c m+ i + m(m + 2)/3 2 c m+2 = 0; 



( m + i)pEH - (AT - m - 2)/3A— 

Cat Cat 



(3) a = 0, a = 67 + p, /3 = 2677, A = 7/9/(2??), p = p — lry: 
Cm - (2TV - 2m + l)Ac m _i - 



(A r_2m + l)p-A^ + ^ 



Cat CAT 

2br]c m+ i + (TV - m)(TV - m + 2)A 2 c m _ 2 
t (TV - m + 1)(TV - 2m + l)Ap + (TV - m)A ( — - A^i 

\CN CAT 

(m + 1)(TV — m)/9 2 — (TV — m)6 7 p + mp— — mXp — 1 

Cat Cat 



Cm— 1 



2bi]pc\ — 2bi]C\ 



c m + 2(m + l)br]pcm+i = 0; 



(4) a = 7 = 0, cc = p, f3 = 2brj, A = 0, p = a: 

CN 



Cm ~\~ 



(TV - 2m - l)a - 



CAT 



Cm - /3c m +i 



(B.4a) 



(B.4b) 



(B.4c) 
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(m + l)(iV — m)a' 2 — ma— + a/3ci — /3ci 

c N 



c m +(m + l)aPc m+1 = 0; 
(5) a = 2am /3 = 2bm p = -67, A = [(2ar? - ^h]/^), /i = 2ar/ - 25 7 : 
c m - (2JV - 2m + l)Ac m _i 



(iV-2m-l)67-A^ + ^ 



(N + 2m + l)(iV - m + l)tr/\ + (N - m)A A 



CN C N 

CN-1 CN 



CN C N 



Cm— I 



+ 



(m + 2)(N - m)b Y + A 



C7V-1 



Cat Cat/ 



+ 2(iV — m)ab r )r\ + 2b 2 jr]Ci + 26r/ci 



c m - 2(m + l)6 2 7?7c m+ i = 0. 
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